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ABSTRACT 

Aims. Our aim is to observationally investigate the cosmic Dark Ages in order to constrain star and structure formation models, as 
well as the chemical evolution in the early Universe. 

Methods. Spectral lines from atoms and molecules in primordial perturbations at high redshifts can give information about the 
conditions in the early universe before and during the formation of the first stars in addition to the epoch of reionisation. The lines 
may arise from moving primordial perturbations before the formation of the first stars (resonant scattering lines), or could be thermal 
absorption or emission lines at lower redshifts. The dilficulties in these searches are that the source redshift and evolutionary state, as 
well as molecular species and transition are unknown, which implies that an observed line can fall within a wide range of frequencies. 
The lines are also expected to be very weak. Observations from space have the advantages of stability and the lack of atmospheric 
features which is important in such observations. We have therefore, as a first step in our searches, used the Odin* satellite to perform 
two sets of spectral line surveys towards several positions. The first survey covered the band 547 - 578 GHz towards two positions, 
and the second one covered the bands 542.0-547.5 GHz and 486.5 -492.0 GHz towards six positions selected to test different sizes 
of the primordial clouds. Two deep searches centred at 543.250 and 543. 100 GHz with 1 GHz bandwidth were also performed towards 
one position. The two lowest rotational transitions of H2 will be redshifted to these frequencies from z~ 20-30, which is the predicted 
epoch of the first star formation. 

Results. No lines are detected at an rms level of 14-90 and 5-35mK for the two surveys, respectively, and 2-7 mK in the deep 
searches with a channel spacing of 1-16 MHz. The broad bandwidth covered allows a wide range of redshifts to be explored for 
a number of atomic and molecular species and transitions. From the theoretical side, our sensitivity analysis show that the largest 
possible amplitudes of the resonant lines are about 1 mK at frequencies < 200 GHz, and a few fiK around 500 - 600 GHz, assuming 
optically thick lines and no beam-dilution. However, if existing, thermal absorption lines have the potential to be orders of magnitude 
stronger than the resonant lines. We make a simple estimation of the sizes and masses of the primordial perturbations at their turn- 
around epochs, which previously has been identified as the most favourable epoch for a detection. This work may be considered as an 
important pilot study for our forthcoming observations with the Herschel Space Observatory. 
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- ISM: molecules - Submillimeter 



1. Introduction 

One of the most important topics in astronomy today con- 
cerns the formation of the first stars and structure formation in 
the Universe. The cosmic microwave background radiation ob- 
served by NAS As Cosmic Background ExploreiQ (COBE) satel- 
lite and the Wilkinson Microwave Anisotropy Probe' (WMAP) 
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* Odin is a Swedish-led satellite project funded jointly by the 
Swedish National Space Board (SNSB), the Canadian Space Agency 
(CS A), the National Technology Agency of Finland (Tekes) and Centre 
National d'Etudes Spatiales (CNES). The Swedish Space Corporation 
was the prime contractor and also is responsible for the satellite opera- 
tion. 

' |http : //lambda . gsf c . nasa . gov/| 



shows small density fluctuations at a redshift of z« 1100 and 
a temperature of ^3000 K. After this epoch of recombination 
the Universe became neutral and entered the cosmic Dark Ages 
which did not end until the formation of the first stars and 
quasars , which are be l ieved to have reionised the Universe at 
z~ 11 jKomatsu et all 120091: lOunklev et al.L l2009l) . During the 
Dark Ages the transition from the small density fluctuations left 
over from the infl ation to completely formed objects took place 
(e.g. lLoebll2008h . The challenging question is how and by what 
means we can reveal this process. 

The main method to obtain information about physical and 
chemical conditions in star forming regions at low redshift is 
by means of atomic and molecular line observations. Molecular 
observations of high redshift objects have also been performed, 
for example bv lBertoldi et alJ (l2003h who found high excitation 
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CO in a Sloan Digital Sky Survey quasar at z-6.4. However, no 
observational evidence exists at all from the cosmic Dark Ages 
at even higher redshifts before or during the earliest epoch of 
the first star formation at predic ted redshifts of ~20-30 (e.g. 
iGloved. l2005t iGlover et all l2008l) . 

Model-dependent theoretical analysis and computer simula- 
tions have therefore been the main tools to follow perturbations 
of different scales in the primordial medium to predict when and 
how the first stars and structure in the Universe evolved. This 
depends on a number of parameters including the chemistry in 
the Early Universe and the properties of dark matter. The chem- 
istry constitutes a very important part in the Early Universe since 
star formation requires a cooling mechanism mainly provided 
by molecules. Dust grains are very important for molecular pro- 
duction and for instance almost all of the molecular hydrogen 
production occurs today on the surface of dust grains. The lack 
of dust and heavy elements in the Early Universe has therefore 
resulted in very low molecular abundances of only a few species 
and therefore difficulties to ex plain the formation of the first stars 
(e.g-lBromm & Larson!, |2003|Gloveil|200l- 

The Standard Big Bang Nucleosynthesis (SBBN) model pre- 
dicts the formatio n of H, D, He and Li a few minutes after the 
Big Bang (e.g. iSteigmanL l2007h . Their primordial abundances 
in this model depend only on the density of baryons which 
has been measured w ith high precision, Qb = 0.0456±0.0015 
(iKomatsu et al. L |2Q09|)- In the expanding and cooling Universe, 
molecular synthesis of for example H2, Hj, HeH^, HD, HD^, 
LiH, and LiH^ could begin as soon as neutral atoms appeared, at- 
taining maximum abun dances around redshifts of 100-400 de- 
pending on species (e.g. lVonlanthen et al.L 12009). The molecular 
primordial abundances depend on many parameters and many at- 
tempts have been made to follow the chemical evolution during 
the Dark Ages with results that often differ by orders of magni- 
tude (e.g. Lem & Shull, 1984; Puv et al., 1993; Galli & Palla, 
1998, 2002; 'LeDDetal., 2002; Black, 2006; Puv & S ignord 
2007; Schleicher et al.. .2008; .Sethi et al.„ .2008; .Signore & Puvl 



2009tlBovino et alil2009h . 



The behaviour of gas in galaxy formation is still an 
open problem, despite the development of simulations, more 
precise semi-analytical models and more complete observa- 
tional data during the recent years. In particular, cooling 
mechanisms in the primordial medium, star formation and 
feedb ack processes ar e far from being completely under- 
stood (ICia rdi & Ferrara. 20 05: Baugh, 2006; Ellis & Silk, 2007; 
lOkamoto ', 2008; De Lucia & Poggianti, 2008). In order to be 
able to discriminate between the different models of structure 
formation it is therefore important to test theory with observa- 
tions from the Dark Ages. 

We do know that during this epoch, the Universe is filled 
with Cosmic Microwave Background (CMB) photons, atoms, 
molecules, ions and electrons. The question is if these species 
have high enough abundances to produce detectable signals. 
Since the average density on large scales is predicted to be 
very low, ne ~ 0.2 cm"^ at z ~ 100 and decreasing towards 
lower redshift, collisions are negligible and excitation is gener- 
ally believed to be dominated by absorption of CMB photons 
followed by spontaneous or stimulated emission in a process 
known as resonant scattering. This process alone is not able to 
produce a signal since CMB photons are nearly isotropic and 
resonant scattering produce an isotropic distribution of photons. 
The only effect of resonant scattering by a primordial medium is 
the damping of CMB primary anisotropics (Maoli et al., 1994). 
However, the gravitational potentials induce peculiar velocities 
of the primordial clouds. Thus, the scattered CMB photons will 



be isotropically distributed in the reference frames of the clouds, 
but not in the CMB reference frame. As a consequence, sec- 
ondary resonant scattering anisotropics will be produced. The 
cross section for the resonant scattering is many orders of magni- 
tudes larger than for Thomson scattering between CMB photons 
and electrons, but it is extremely frequency dependent, peak- 
ing at the rotational and rotational-vibrational frequencies of the 
scattering species. Resonant lines will therefore be produced by 
these moving primordial clouds ( Ma oli et"al],ll996h . 

Since the original idea bv iDubrovichI d 1 9771) . inspired by ear- 
lier work on CMB fluctuations by Zeldovich, a number of papers 
have alre ady analysed the possibility of detecting the resonant 
lines (e.g. iMaoli etaLl Fl 994. 1 996; Dubrovich et al.. 2008) and 
two pre vious observational attempts have been made with out 
success ( de Bernardis et al.Lll993l : lGosachinskii et am2002h . 

However, possible signals from the Dark Ages may not only 
arise from resonant scattering since the growing potentials may 
further induce adiabatic heating as well as shock heating of col- 
lapsing primord i al per t urbations (e.g . Barkana & Loeb, 2001 ; 
'Santoro & Shul l2006t iJohnson & Bromm. 2006; Grei fetaT , 
2008b). This produces thermal emission or absorption from 
denser than average or high temperature regions. The first stars 
are predicted to be very massive >1OOM0 and formed either in 
isolation or at most as a small stellar multi ple in minihalos (e.g . 
lAbel et al.Ll2000l;lBromm & Larsonl l2004; Bro mm et"alll2009l) . 
These hot stars must have emitted enormous amounts of ener- 
getic radiation. By the time of reionisation, large HII regions are 
predicted to have g rown around clusters of these hot and bright 
stars (e.g. lGreif et a l.. 2008a). This radiation will first destroy the 
molecules, but in the cooling ejecta of the primordial supernovae 
molecules will later be efficiently produced (Cherchneff & Lilfvl 
2008|) for instance H2 due to the high degree of ionisation. 
However, many considerations have to be taken into account to 
predict all these types of emission/absorption, hence leading to 
uncertain results. 

We emphasize that every possibility to obtain observations 
from the Dark Ages should be seriously considered. A detection 
would be extremely important and would give direct evidence of 
a very high redshift protostructure within a wide range of angu- 
lar sizes, from arcsecond to arcminute scales, depending on the 
cosmological model. These small scales are not observable in the 
CMB anisotropics due to the finite depth of the last scattering 
surface. The horizon at the last scattering surface corresponds 
to 240 Mpc today. The PlanclQ satellite, launched on May 14, 
2009, will measure the CMB power spectrum for angular scales 
greater than 10' (40 Mpc). Smaller scales, at cluster and galaxy 
sizes are unobservable due to photon diffusion damping asso- 
ciated with the non-instantaneous nature of the recombination 
process at the last scattering surface (Silk, J_968,). 

Any primordial spectral line will also provide a unique test 
of nuclear synthesis at high redshift, give new information about 
the reionisation of the universe, probe the chemistry throughout 
a wide range of redshifts, the heating and cooling processes as 
well as the dynamics of the primordial clouds before and dur- 
ing the gravitational collapse of a protostructure. Even a non- 
detection could give very valuable information in order to con- 
strain all these issues if the noise level is low. 

Given the expected weakness of the lines and our ignorance 
of the frequencies at which the transitions will fall, ground-based 
observations are complicated by the terrestrial atmospheric lines. 
The ozone molecule, O3, including all the isotopes '^O, '^O and 
'**0, emits numerous lines. Moreover, most of the vibrationally 
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excited lines can not be foun d in catalogues. Observati ons with 
the IRAM 30-m telescope (Ide Bernardis et al.L Il993h experi- 
enced difficulties from these lines at a low level. Interferometer 
observations are therefore interesting for searches like this, due 
to the many advantages for rejecting both broad and narrow band 
atmospheric emissions. Ground based single dish observations 
can, however, also provide valuable information especially at 
low frequencies, even though the foreground and atmospheric 
radiation poses a large problem. Another possibility is observa- 
tions from space, which have the advantages of stability and lack 
of atmospheric features. 

As a first test, we have therefore chosen to use a satellite to 
search for and put upper limits on primordial signals from the 
Dark Ages. The satellite is required to have tunable receivers in 
order to cover a broad spectral band, essential for the exploration 
of a wide redshift range. In this work, we have used the only 
available satellite at the tim e with the above mentioned require- 
ments: the Odin satellite (.Nordh et al.L l2003l) . This satellite has 
the unique capability to cover a broad spectral band with a high 
channel spacing of about 1 MHz using tunable single sideband 
(SSB) receivers. This has allowed us for the first time to search 
for redshifted molecular hydrogen from the predicted time of the 
formation of the first stars or earlier. 



2. The Odin observations 

The main problems in the searches of primordial lines are the un- 
known but expected weak amplitudes from unknown species in 
unknown sources at unknown redshifts. This implies that the fre- 
quency of the lines can be anywhere in a wide frequency range. 
Thus we need to cover a large frequency band to probe a wide 
range of redshifts for many possible species and transitions. The 
sources should on the other hand be found in every direction on 
the sky even though they most likely have a clumpy distribution 
and also are located at diff'erent redshifts. This implies that we 
do not know the evolutionary state, size or density of the pertur- 
bations. 



2. 1. Interesting species 

Our approach in all observations has been to perform spectral 
line surveys in order to cover a wide redshift range. This is nec- 
essary since our objects are located at an unknown redshift and 
the expansion of the Universe causes the frequency of an emitted 
photon to be redshifted as 



Vobs 



1 +z 



(1) 



where vq is the rest frequency, Vobs the observed frequency of 
the transition, and z is the redshift. When we choose an observa- 
tion frequency, we thus determine the emitting redshift for each 
transition given by Eq. ([1]). Each transition may also be detected 
at different frequencies, emitted from several objects at differ- 
ent redshifts along the line of sight. The strong frequency de- 
pendence of the scattering and emission processes gives rise to 
lines superposed with the black body continuum spectrum of the 
CMB. Each line is associated with a moving primordial pertur- 
bation, exactly in the same way the Lyman a absorption lines are 
associated with neutral hydrogen clouds which absorb the light 
from bright quasars and galaxies. 

We search for any atom or molecule that may be present 
in the high redshift Universe. This includes mol ecules that 
form from the primordial elements H, D, He and Li (ISteigmanL 



l2007h such as: H2, H+, H2D+, HD, HD+, HeH+, LiH and 
LiH^. There are also suggestions that lines from neutral hy- 
drogen H and helium He could be seen during the time of re 



spect i ve recombination and at later times ( Rubino-Mart m et al. 
'2006UChluba & Sunvaevil2006t lSethi et al L2007-,Chluba et al. 
2007; Rubino-Martm etal., 200^ Prior to the reionisation at 
z > 6 the universe was, however, opaque to for instance Lya 
radiation! Becker et al., 2001). 

Previous searches for primordial resonant lines were mainly 
focused on LiH due to its high dipo le moment and a possibly 
high abundance dde Bernardis et all [1990). New ab initio cal- 
culations of reaction rates of LiH destruction in the early uni- 
verse show, however, that t he LiH abundance is too low to allow 
a detection jBovino et al. I l2009t) . Another interesting molecule 
with a high dipole moment is HeH^ which, together with a much 
higher abundance than LiH, increases the possibility of a detec- 
tion. 

Moreover, non-standard BBN models predict heavy nuclei 
such as carbon, oxygen, nitrogen, and flourine, which may 
have created molecules like CH , OH, NH, HP, and their re- 
spective molecular ions (jPuv et al.. ,2007: Carnpos et al., 20071 
IVonlanthen & Puvi 120081: IVonlanthen et al.L l2009h . For exam- 
ple the abundance of the CH molecule is found to be higher 
than the H2D or HD abundances in the inhomogeneous non- 
standard BBN models. In eff'ect, a number of possible mecha- 
nisms have been suggested to generate density perturbations in 
the early universe which could survive until the onset of primor- 
dial nucleosynthesis. Such inhomogeneities in the baryon num- 
ber may result from non-equilibrium processes in the Big Bang, 
for example, occurring during a putative first-order quantum- 
chronodynamics (QCD) phase transition at ~ 100 MeV or during 
the electroweak symmetry breaking at ~100GeV. Assuming the 
existence of such surviving density fluctuations, some authors 
have shown that - in the framework of these non-standard BBN 
scenarios - the altered primordial nucl eosynthesis could lead 



to the sy nthesi s of hea v y elements (e.g.lApp 



egateet all 11987 



Rauscher et all 119941: iKurki-Suonio et all 119971: iLara et" 
12006 ). Note also that these models satisfy the observed primor- 
dial elemental light abundances. 

In addition, regardless of BBN model, the first genera- 
tion stars of a primordial composition formed at z ^ 20-30 
might quickly have polluted the medium with small amounts 
of metals with subsequent che mistry (e.g. Mackev et alj 



' 20031 



equent chi 

jBrornm & Larsonl l2004t ICherchne ff & Dwek. ' 200', , 

Sm ith et al 11009). Several species have already been observed 
at high redshifts. In additi on to the high excitation CO ob- 
servations at z - 6.4 ( Bert oldi et al.L l2003h this molecule has 
also been obs erved in for example a damped Lyman-a system 
at z = 2.4 (Srianand et al.L I2b08i) and in a field containing an 
ove r-density of Lyman break galaxies at z = 5. 1 (IStanwav et al.L 
2008). A mas sive CO reser voir has also been detected at z = 3.9 
(Papadopoulo s et all 1200 ih . An un usually high amount of neu- 
tral hydrogen was found bv iFrve et al.l (|2008) in a 14Mpc re- 
gion surrounding a young galaxy at z =4.9. A strong detection 
of molecular absorption bands, including H2 and CO, was also 
recently observed from gas within the host galaxy of the gamma- 
ray burst 080607 at z = 3.0363 (Prochaska et al., 2009). Several 
absorption features are still unidentified, but a H2 column density 



of 10 



,2I.2±{).2 



cm was inferred with an excitation temperature 



of 10-300 K. The highest redshift object o bserved today is th e 
gamma-ray burst of 23 April 2009 at z = 8.2 (iTanvir et al.Ll2009l) . 
The long duration gamma-ray bursts are believed to be primor- 
dial supernovaes originating from t he super-massive first gen er- 
ation stars, i.e. population III stars dWooslev & BloornLl2006h . 
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The continuum of has also been suggested to give rise 
to a detectable signal in terms of a decrease in the CMB spec- 
trum which may superimpose a detectable absorption feature on 
the CMB. This could be detectable with the Planck satellite, al- 
though the strength of this e ffect has been estimate d to quite dif- 
ferent values (Black, 2006; Schleicher et al.Ll2008h . 

Besides producing molecular lines, resonant scattering has 
an additional interesting feature that is regardless of peculiar ve- 
locities, this process has the potential of reducing the power of 
the CMB primary anisotropies i f the molecular abundances are 
high enough (Ma olietal.Lll99l . Planck may be able to detect 
also this decrease in power There are also several suggestions 
to use differential measurements of the broad band CMB an- 
gular power spectru m observed with Planck in the sea r ch for 
resonant lines (e.g. | Basu et al. I l2004t iDubrovichetaTl llool 
[Schleich er et al.L 12008 ). To reach the required sensitivity of the 
order of yuK, the main problem to extract these weak signals with 
this approach, will be the galactic and extra-galactic foreground 
emissions. 

Since the most abundant species is neutral hydrogen, a very 
interesting and promising transition is the 21 cm line. Such diffi- 
cult observations must, however, await the future Low-frequency 
ArrayPi (LOFAR) and the proposed Square Kilometer Array 
(SKAj3 which will be able to detect the redshifted 21 -cm HI 
transition between z~6-ll or at even higher redshifts (e.g. 
[Hogan & Rees, 1979; Furlanetto et al., 2009). Both instruments, 
however, face severe difficulties when subtracting the astrophys- 
ical foreground contamination wh ich is several orders of magni- 
tudes stronger than the 21cm line dBowman et al.Ll2009l) . 

Most other atomic transitions than the 21 cm line have much 
higher required excitation temperatures than the molecular tran- 
sitions and, except for H, He and Li, must in addition have been 
produced either by the first stars or in a non-standard BBN. 
They are therefore only expected to be seen at redshifts dur- 
ing the epoch of first star formation or some time after. As a 
first step, we have therefore concentrated our searches on opti- 
cally thin primordial molecular lines. The redshifted frequencies 
from z ~20-30 of the lowest transitions of ortho- and para-H2 
fall at Odin frequencies around 500 GHz. These lines could be 
observable in absorption towards the CMB or more likely to- 
wards for example hot HII regions. The detectability of H2 ro- 
tational lines in emission associated wit h the formation of the 
first stars has previously been discussed (lOmukai & KitavamaL 
[2003; Kamaya & Silk, 2003|^ in addition to the ro- vibrational H? 
emission lines (Mizus awa et al.L [20041) . In addition, a number of 
lines of for example Hj, HiD^, HD, HD^, and HeH^ can also 
be seen from this epoch. 

At redshifts higher than about z> 150-200 the kinetic matter 
temperature is expected to be the same as the radiation tempera- 
ture of the CMB due to Compton scattering, evolving as 



T'cMB - T]^ - Tod + z) , 



(2) 



where Tq = 2.725 K is the cosmic microwave background radi- 
ation (CMB) today, as m easured by the FIR AS instrument on- 
board the COBE satellite (iFixsen et alill996h . During this epoch 
only resonant scattering may occur since temperature differences 
are required to produce thermal spectral lines (see Eq. l23T l. At 
redshifts below 200, the matter temperature evolves following 
an adiabatic expansion which i mplies a kinet ic temperature de- 
crease faster than for the CMB dGloved [20051) 



(3) 



These temperature differences will force the excitation further 
away from equilibrium with the radiation, and should enhance 
spectral absorption distortions at redshifts < 150-200. The gas 
will eventually be heated by the accretion and collapse phases of 
the growing perturbations, and eventually by the first stars and 
quasars which will reionise the whole universe at z> 6. This will 
enhance spectral emission features. 

The coupling of the excitation temperature to the radiation 
or matter temperature differs with species. Polar molecules like 
HeH^ probably continue to have Tex coupled to the CMB tem- 
perature, while H2 may remain more strongly coupled to the 
matter temperature. There are also situations where the excita- 
tion temperature may depend on molecular formation processes, 
by which for instance primordial H2 may remain superthermal 
relative to the matter temperature. 

The strength of a spectral line depends, in addition to the 
number density and dipole moment, on the ratio of the upper 
state energy (in Kelvin) and excitation temperature (see Ea.[T3ll. 
Rotational-vibrational Unes, with upper state energies of thou- 
sands of Kelvin, are therefore only expected for very high red- 
shift sources with a high temperature or from collapsing high- 
temperature regions at lower redshifts. In such regions the den- 
sity may be increased above average, and the increased collision 
rate in these regions will drive the level populations and the ex- 
citation temperature towards the high matter temperature. 

A detection of several spectral lines is required to secure a 
molecular identification and a redshift of an emitting source. 
The spectral distance between rotational lines is, however, rather 
wide, even though it becomes increasingly more narrow for 
highly redshifted sources observed below 100 -200 GHz. It is 
thus more likely that a single transition from each object would 
be seen in the Odin observations. The rotational-vibrational tran- 
sitions would be more closely spaced, but have a lower probabil- 
ity of detection. 

2.2. Observational strategy 

Neither the amplitude nor the line widths are known and there- 
fore we use high spectral resolution to enable detection of nar- 
row lines, and later re-bin the data to lower resolution. The 
choice of spectral resolution affects the root-mean-square noise 
level Trms as can be seen from the radiometer formula 



K 



(4) 



- sys 



^ |http : //wm . lof ar . org/] 

* [http : //www ■ skatelescope . org/| 



where A/ is the frequency resolution of the spectrometer, t is the 
on-source integration time, T,,., is the noise temperature from 
the whole system, which includes the receiver noise, and A' is a 
constant which depends on the observation strategy (K- V2 for 
a switched receiver). 

The expected signals most likely originate from every direc- 
tion in the sky, even though they will be stronger from rare high 
density peaks. However, we needed to observe positions with 
as low contamination as possible from any known source, and 
in addition they had to be observable with the Odin satellite for 
long periods. For these reasons we have observed towards two of 
the WMAP hot spots in the CMB radiation, out of the Galactic 
plane. The coordinates are listed in Table [1] 

Two different sets of observations were performed with dif- 
ferent observation strategies as described below. The first survey 
took place in 2004 and the second one in 2006-2007. In addi- 
tion, we performed two sets of deep searches at one frequency 
setting towards one position in 2007 and 2009. 



CM. Persson et al.: The first spectral line surveys searching for signals from the Dark Ages 



5 



Table 1. Coordinates of observations. 



Position 


R.A. (J2000) 


DEC. (72000) 


Hot Spot 1 


05"26'"00'0 


-48° 30: Off.'O 


Hot Spot 2 


05''09'"36'0 


-43° 24: oa'o 



The Odin 1.1m offset Gregorian telescope ha s a circular 
beani at 557 GHz and a FWHM beam width of I'A (iFrisk et al.L 
l2003h . This corresponds for example to 540 and 78 kpc at z = 10 
and z- 100, respectively (see also Table IC2l on-line material). 
Being outside the atmosphere, and with an exceptionally high 
main beam efficiency, Ty^b = 0.9, our intensity calibration is very 
accurate. The calibration procedure (the chopper wheel method) 
of the Odin satellite is described in Olberg et a n (l2003h . The in- 
tensity scale in the figures is expressed in terms of antenna tem- 
perature T^. The reconstructed pointing offset was <15" during 
most of the time. 

Two different single sideband (SSB) receivers, with a typi- 
cal side-band suppression of >20dB, were used simultaneously 
in combination with one auto-correlator (AC) and one acousto- 
optical spectrometer (AOS). The bandwidth is 1.040 GHz for 
the AOS, and 690 MHz for the AC, with a channel spacing of 
0.62 and 1.0 MHz for the AOS and AC, respectively, and thus 
Ay/ V « 2x10"^ at V = 545 GHz. The Odin average system temper- 
ature is around 3 300 K (SSB) and therefore many hours of ob- 
servations are demanded to produce sensitive observations. The 
resulting noise levels may very well be too high to allow a detec- 
tion, but the intensity of the lines is really not known and we have 
therefore taken a pure observational approach. If no detection is 
obtained, the resulting upper limits will be used as input to the 
noise levels required in future observations. This approach also 
applies to the chosen frequency resolution which is connected to 
the noise level and final covered bandwidth. 

In each Local Oscillator (LO) setting, the resulting AOS 
spectra have a very stable baseline except at the band edges 
where the calibration spectra exhibit a steep increase. The noisy 
edges are excluded, and then a first or second order polynomial 
baseline is subtracted before we join all spectra together to form 
a contiguous spectrum. The AC consists of 7 bands of 100 MHz 
each in each setting. We subtract a linear baseline in each band 
before we align and average all spectra. All averages are then 
joined. 



2.2.2. Observations during 2006/07 

In the next attempt during nine weekends in winter 2006/07 and 
one weekend in Aug 2007 (in total 464 orbits) we changed our 
observational strategy. We wanted to lower the noise compared 
to what we obtained in the 2004 observations, but still needed 
to cover a wide frequency band. Since the sizes of the primor- 
dial objects are unknown, this time we also wanted to test dif- 
ferent spacings between the signal and reference positions in 
order to avoid subtraction of the signal if present in both po- 
sitions. A compromise of the above requirements resulted in the 
scheme shown in Fig.[T] Four positions, A-D, were observed 
where position A is the previously Hot Spot 1 observed during 

5' 15' 10' 



• • • • 

A B CD 

Fig. 1. Observation strategy during the second spectral survey 
2006/07 (also described in Table |2|. The Odin beam has a 
FWHM width of 2:1. 

2004. The minimum offset was determined by the Odin beam 
size of 2'.l, and thus we chose the following offsets: -5', -15', 
and -10' in DEC, keeping the same R.A. During 10 weekends 
of observations (~40 orbits/weekend) we observed all four posi- 
tions with one frequency setting each weekend. Table|2]describes 
how we performed the observations: 10 orbits towards each pair 
of positions were observed in a sequence. First towards position 
A with B as reference (AB), then B with C as a reference (BC), 
C with D as reference (CD), and finally D with A as a reference 
(DA). Thus, in total we have ten orbits/setting/position for ten 
different LO settings. 

Table 2. Observation strategy 2006/07: every weekend we ob- 
served four positions, A-D, in a sequence using ~10 hours to- 
wards each position with one frequency setting (also shown in 
Fig.D. 



Sig Ref Spacing 



A 


B 


5' 


B 


C 


15' 


C 


D 


10' 


D 


A 


30' 



2.2.1 . Observations during summer 2004 

Our first goal was to cover as wide a frequency band as pos- 
sible. We were awarded 337 orbits (1 orbit~l hour of obser- 
vation) and performed independent position switching observa- 
tions towards the two positions listed in Table [1] using refer- 
ence positions offset by -45 arcminutes in declination. Using si- 
multaneous observations of the AOS and AC, 15.5 GHz were 
covered in steps with each receiver: AC ~563-578 GHz and 
AOS ~547-563 GHz, giving at total observed frequency band 
of 31 GHz covering the full 547 -578 GHz band. We spent 5 
orbits for each LO setting and the step size in frequency was 
0.5 GHz. The settings and steps were the same as in the first 
spectral line surve y performed by the Odin satellite towards t he 
Orion KL nebulae dOlofsson et al.Ll2007HPersson et al.Ll2007h . 



As in the first observation run in 2004, we performed si- 
multaneous observations with the AC and the AOS, but cover- 
ing other frequency ranges. The final results are spectral sur- 
veys in the frequency ranges 542.0-547.5 GHz (AOS) and 
486.5 -492.0 GHz (AC) towards four positions. 

2.3. Deep searches at 543 GHz during 2007 & 2009 

In addition to the spectral surveys we also wanted to perform a 
deep search to lower the noise even further. In 2007 we there- 
fore used 40 orbits towards HotSpotl position A with B as ref- 
erence position with 543.250GHz as centre frequency with the 
AOS and 490.250GHz with the AC. This produced an rms level 
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of about lOmK over the 1 GHz band with a channel spacing 
of 1 MHz, and a possible 4cr detection at 543.1 GHz. To val- 
idate the possible detection, we performed additional observa- 
tions in April 2009 using 96 orbits and then changed the centre 
frequency of the AOS to be 543. 1 GHz. 

3. Results 

No lines were detected and we are thus limited to set upper lim- 
its of possible signals. We choose a 5cr noise level as an upper 
limit since the peak noise level is always lower than this value 
if the noise follows a Gaussian distribution and has independent 
velocity channels. In this case, the probability that an observed 
5cr line would be a noise feature is less than 10"^. 

The noise is rather stable over the covered band of all sur- 
veys, with a few exceptions where we had loss of observation 
time. We therefore measure the Icr noise level across the to- 
tal band covered, towards each position. We have re-binned the 
AOS data to the AC channel spacing of 1 MHz, and then re- 
binned both AOS and AC to channel spacings of 4 and 16 MHz, 
with an example shown in Fig.|2] This rebinning should lower 
the noise by a factor of 2 and 4, respectively, if the noise 
is Gaussian distributed. However, our measurements do suffer 
from baselines at a low level which prevents the noise to fol- 
low Eq. (|4]i exactly. All Icr results, except the deep searches, are 
found in Tables|3]and|5] which also include the peak intensity of 
the noise given in terms of the measured cr. 



T 



n 

MHz 



Note, that we do not know if the signals we seek are absorp- 
tion or emission lines, or if they are located in the signal or the 
reference position. 

3.1. Comments on results from the 2004 observations 

The observations from summer 2004 resulted in a wide fre- 
quency band, but with a rather high noise level. In Table [3] we 
see that a typical Icr is about 38(85), 23(45), and 15(26)mK 
for the AOS(AC) observations with 1, 4 and 16 MHz channel 
spacing, respectively. 



Table 3. Upper limits from observations during 2004 using 1, 
4 and 16 MHz channel spacings. The limits are for the data ob- 
served with the AOS(AC) and for frequencies between 547-563 
(563 -578) GHz. 



1 MHz channel 

Position Icr 
[mK] 


spacing 
Peak Intensity 
[0-] 


HotSpot 1 37 (88) 
HotSpot 2 40 (83) 


3.8(4.0) 
3.6(3.9) 


4 MHz channel 


spacing 


HotSpot 1 22 (46) 
HotSpot 2 24 (44) 


3.1(3.3) 
3.2(3.4) 


16 MHz channel 


spacing 


HotSpot 1 14(27) 
HotSpot 2 13(25) 


2.9(3.1) 
3.0(3.2) 



-t-> 

• I— t 

M 

ti 

CD 
-i-> 

iB 0.1 

0.05 



-0.05 



-0.1 















H \- 

1 


h 


4 MHz : 

1 . = 


16 MHz 



542 



543 



544 



545 



546 



547 



Frequency [GHz] 

Fig. 2. Examples of the resulting spectra: 2006/07 observations 
towards Hot Spot 1 (position A with position B used as refer- 
ence) is shown with three different channel spacings: 1, 4, and 
16MHz with corresponding Icr noise levels of 20, 13, 7 mK. 



Figures lA.ll and IA.2I (on-line material) show the original 
AOS spectra used to produce Fig.|2l The intensity is plotted vs. 
channel numbers and shown for all ten frequency settings. The 
intensity scale in all figures are expressed in terms of observed 
antenna temperature Ta*. 



3.2. Comments on results from the 2006/07 observations 

Due to the new scheme the analysis of these observations is 
somewhat different from that of the 2004 observations. As seen 
from Fig.[I] each pair of observations is sensitive to different an- 
gular scales of the perturbations. As an example, if we assume 
that the signal is in position A and use B as reference, this com- 
bination is sensitive to a perturbation size between 2-5'. If the 
primordial cloud is larger than 5' the signal will be present in 
both A and B and will disappear in the resulting AB average. 
We are, however, also sensitive to larger sizes of the perturba- 
tions in A since we have observations towards position D with 
A as a reference. If we switch and use A as signal and D as ref- 
erence, we will have an AD average which is sensitive to scales 
between 2-30'. Therefore, we can average both AB and AD 
together to further lower the noise with a sensitivity of scales 
2-5'. For larger scales we can only use AD. A summary of all 
average combinations for each position with corresponding an- 
gular scale sensitivity is given in TableU] 

This observation strategy is also a way to check if the signal 
is in the signal or reference position. For example, if there is a 
signal in position A, and the size of the perturbation is 2 - 5', the 
signal-to-noise (S/N) will increase by V2 if we average AB and 
AD. If the signal is in B, the signal will be lowered by a factor 
of two in the AB-l-AD average. A signal in the B position would 
also be seen in the BC average with V2 increase of the S/N in 
the BA+BC average. 

In Table |5] we note that the noise level is considerably lower 
than in the 2004 observations. 
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Table 4. Observational strategy 2006/07: different combinations 
of signal and reference positions are sensitive to different angular 
sizes of the primordial perturbations. 



Position 




Angular sensitivity 6 






2-5' 


5-10' 


10-15' 


15-30' 


A 


AB+AD 


AD 


AD 


AD 


B 


BC+BA 


BC 


BC 




C 


CD+CB 


CD+CB 


CB 




D 


DA+DC 


DA+DC 


DA 


DA 



3.3. Comments on results of the deep searches 

The noise level in both deep searches are considerably lower 
than the 2004 observations and also lower compared to the sur- 
veys in 2006/07. The deep search in 2007 resulted in a noise 
level of lOmK with 1 MHz channel spacing. Rebinning to 4 
and 16 MHz lowers the noise to 6.5, and 3.5 mK, respectively. 
The noise level in the 2009 deep search produced an even lower 
noise level of 6.7 mK with 1 MHz channel spacing, and 4.0, and 
2.1 mK when rebinning to 4 and 16 MHz, respectively. 

The possible detection in the 2007 deep search did unfor- 
tunately not show up again in the 2009 deep search as seen in 
FigH The origin of this feature is difficult to explain but may 
arise from unknown technical issues with the satellite or the re- 
ceivers. 



■-0.05 
0) 



-0.05 



' 1 ' 1 
. Deep Searcli 2007 

1 1 1 — 


' 1 ' 1 

4 MHz 
^ 1 1 [ 


. Deep Search 2009 

1 , 1 


4 MHz 

1 , 1 



543 



543.2 



543.4 



543.6 



Frequency [GHz] 

Fig. 3. Deep searches from 2007 and 2009 observations towards 
the Hot Spot 1 (position A with position B used as reference) is 
shown with a channel spacing of 4 MHz with Icr noise levels of 
6.5 and 4.0 mK, respectively. 



4. Analysis 

In order to interpret our observations in a cosmological con- 
text we begin with a summary of the cosmological background 
and the tools by which we estimate the amplitude of absorp- 
tion and resonant lines in Sect. |5] as well as their line widths. 
This is model dependent and we adopt the currently favoured 



hierarchical ACDM cosmology with a cosmological constant in 
which cold dark matter dominates the evolution of structu re. The 
latest results from five-years of WMAP data ( Komatsu et all 
l2009h are used, with a Hubble constant Ho = 70km s ' Mpc 
a dimensionless matter density parameter 0.274, a dark 
energy density parameter Qa - 0.726, and a baryon density 
Qb = 0.0456. This implies a flat geometry where Qm + = 1- 

The proper length L(z) of a primordial perturbation can be 
estimated by using the small-angle limit of the relation between 
the angular size 6 that a cloud subtends on the sky and the 
angular-diameter distance DpXz): 



L(z) 



[rad]. 



Da{z) 

In a flat Universe, this distance is described by 

dz 



£>a(z) 



(1 



+ z')HoL 



VnM(l+z)3+n^ 



[Mpc] 



(5) 



(6) 



where c is the speed of light. Assuming a spherical geometry, the 
proper length L(z) can be used to calculate the mass of a density 
perturbation Ap/p following the expansion of the Universe 

M = y ^ Pci-( 1 + 1 + Ap/p) , (7) 

where pe,- = l.SSxlO"-''/!^ gem ^ is the critical density of the 
Universe at present time with /i = //o/100, and where Ap/p 
is initially very small but increases with decreasing redshift. 
Equations (|5]l - Q then give a relation between angular size and 
mass. 

The standard structure formation model predicts the forma- 
tion of gravitationally bound systems from tiny density perturba- 
tions via gravitational collapse. The hierarchical model predicts 
that the smaller perturbations formed first and then merged or 
accreted gas to form even more massive objects. The WMAP5 
results are consistent with an ep och of reionisation at Zra = 
10.9 +1.4 (iKomatsu et al.L l2009h . Later reionisation epochs at 
z~6, are suggested by other kinds of experiments (e.g. lFan et al.L 
2006). This implies that the first stars must have formed at z>6. 
They are predicted to form in over-dense dark matter regions o f 
10^-10^ Mo at redshifts of about 30-40 (e.g. iGloveil l2005h . 
These mini-halos may provide a significant, if not dominant , 
contribution to reionisation at lower redshifts (Chou dhurv et al.L 
2008). Mass perturbations >1O'^M0 are predicted not yet to 
have reached their collapse phase. 

In the linear regime, the density contrast Ap/p is predicted to 
grow linearly with the scale factor as 1/(1 +z). Sufficiently dense 
perturbations reach a turn-around point at which their gravity 
counterbalances the expansion. Then they enter a collapse phase 
during which the density and temperature quickly increase. 

Note that during the linear phase, the only predicted possi- 
ble signal is by the resonant scattering process, while emission 
and absorption lines will arise in collapsing perturbations or with 
high temperature regions for example from the first supernovae 
and their subsequent hot HII regions as background radiation. 
The amplitude and line width of all lines depend among many 
variables on the dynamics of the primordial clouds, which con- 
veniently can be divided into the above mentioned phases: the 
linear phase, the turn-around phase and the collapse phase. 



4.1. Linear phase: 

In the linear phase, the line width will depend on the proper 
length L{z) of the object. This size will occupy a redshift interval 
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Table 5. Upper limits in mK as a function of angular size on the sky from the spectral surveys observations during 2006/07 (deep 
searches are not included) using 1, 4 and 16 MHz channel spacings. The limits are for the data taken with the AOS(AC) and for 
frequency ranges 542.0-547.5 (486.5 -492) GHz. 



1 MHz channel spacing 



Position 


2-5' 


Peak intensity 


5-10' 


Peak intensity 


10-15' 


Peak intensity 


15-30' 


Peak intensity 




lo- [mK] 




Icr [mK] 


[crl 

J 


Icr [mK] 




Icr [mK] 


[crl 

J 


A 


16 (23) 


3.0(3.3) 


21(31) 


3.2(3.4) 


21 (31) 


3.2(3.4) 


21(31) 


3.2(3.4) 


B 


17 (26) 


3.1(2.9) 


23 (34) 


3.3(3.2) 


23 (34) 


3.3(3.2) 






C 


19(26) 


3.3(3.0) 


19(26) 


3.3(3.0) 


25 (35) 


3.1(3.2) 


- 


- 


D 


15 (25) 


3.1(2.9) 


15(25) 


3.1(2.9) 


20(33) 


3.2(3.1) 


20(33) 


3.2(3.1) 










4 MHz channel 


spacing 








A 


10(15) 


3.0(2.8) 


13(20) 


3.0(2.9) 


13(20) 


3.0(2.9) 


13(20) 


3.0(2.9) 


B 


11(17) 


2.7(3.0) 


15(23) 


2.9(3.1) 


15 (23) 


2.9(3.1) 






C 


12(16) 


3.0(2.9) 


12(16) 


3.0(2.9) 


16(22) 


2.9(3.0) 






D 


10(15) 


2.9(3.1) 


10(15) 


2.9(3.1) 


13(20) 


3.1(3.0) 


13(20) 


3.1(3.0) 










16 MHz channel 


spacing 








A 


5(8) 


2.7(2.8) 


6(11) 


3.0(2.9) 


6(11) 


3.0(2.9) 


6(11) 


3.0(2.9) 


B 


6(9) 


2.8(3.0) 


8(12) 


2.7 (3.1) 


8(12) 


2.7(3.1) 






C 


6(11) 


2.7(2.8) 


6(11) 


2.7(2.8) 


8(14) 


3.0(2.9) 






D 


5(8) 


2.9(3.0) 


5(8) 


2.9(3.0) 


7(11) 


2.8(3.1) 


7(11) 


2.8(3.1) 



Az in the Hubble flow, and, assuming a spherical geometry and 
that every part of the cloud moves w ith the same pecu liar veloc- 
ity, we have (e.g. lMaoli et al.Lll996tlDubrovich et al.Ll2008h 



Av 



Az HoLjz) I 

- — = v^^m(1 + zy + 

V I + z c 



(8) 



where Av is the Doppler line width. The line width is thus de- 
pendent on the size of the perturbation. 

A relationship between line width and mass in the linear 
phase can now be found using Eq. Q with Ap/p« 1 and Eq. 
The line widths increase with redshift and mass and are very 
broad, Av/v~10"' - 10"^, implying line widths of the order of 
a few thousand kms ' (Table |C3] on-line material). Note, that 
this is only true assuming that the density contrast is very small. 
When the perturbations grow the line widths will start to devi- 
ate from Eq. (O and become increasingly more narrow. Another 
addition of uncertainty to the line width is our homogeneity as- 
sumption of no substructure within the cloud. Within each mass 
at its turn-around, there will be smaller mass perturbations which 
already have reached their respective turn-around. The signals 
from these regions will be very naiTow and superposed on the 
broader lines from the non-collapsed regions. As a first order 
approximation in our analysis, we will therefore use a line width 
of 500 kms ' in our predictions of HeH^ resonant lines in the 
linear phase in Sect.|5] 

The intensity of the resonant lines is est i mated to be (e.g . 
de Bemardis et all 119931: iMaoli et all 119961: iDubrovich et all 
200i 



A/ _^ Vp 
- (l - e ) (3 - ay) — cos ( 

^CMB c 



(9) 



where A/ is the observed intensity, /cmb the CMB intensity, ay 
is the spectral index, Vp is the peculiar velocity of the perturba- 
tion, cos 6 is the cosine of the angle with respect to the line of 
sight of the peculiar velocity, and r is the opacity of the transi- 
tion as calculated in the rest frame of the cloud. The determining 
factor resulting in a positive or negative sign of A/// are the 



cos 6 factor The only way to produce resonant scattering "emis- 
sion" lines is therefore to have a source that moves along the line 
of sight towards us, while resonant scattering "absorption" lines 
will appear when the source moves away from us. Note, that the 
resonant lines do not require temperature differences, but only a 
moving cloud, CMB photons and the scattering species. 

The spectr a l inde x of the CMB ay = (v//)(d//dv) is derived 
in lMaoli et alj (Il996h where ly - By(TcMB) and has the general 
expression 



h v/kT 

I - g-hv/kT 



(10) 



The peculiar velocity describes the motion of primordial per- 
turbations with respect to the Hubble flow. In the ACDM model 
these perturbations evolve due to potential gradients and increase 
with time as ( Longair. 2008) 



(11) 



The most commonly used peculiar velocity at present time 
for the size of a galaxy cluster is Upo ~ 600 kms"', derived 
from the CMB dipole anisotropy, which together with our spa- 
tially flat ACDM model, is considered to reproduce the char- 
acteristics of the large-s cale matter distribution at low red- 
shifts (e.g. iJenkins et al.L [T998). Smaller mass perturbations 
most likely move with a factor of 2-5 higher peculiar velocity 
(IDubrovich et al.L[2008h . The present ratio of the peculiar veloc- 
ity to the speed of light is ~2x 1 0"^ and decreases with increasing 
redshift. This is a very limiting factor in the search for resonant 
lines. 

The scattering efliciency of resonant lines is very frequency 
dependent and depends mainly on two parameters: the number 
density of the species and the spontaneous transition probability 
described by the Einstein coefficient A„/, where the subscript ul 
refers to the upper and lower levels of the scattering species. 
Those species which have a high product of these quantities are 
of special interest. 



CM. Persson et al.: The first spectral line surveys searching for signals from the Dark Ages 



9 



The density of the species depends on redshift and can be 
described in terms of the total density of hydrogen and the frac- 
tional abundance of species x as - n jnu. For simplicity, as a 
first approximation we assume that the density is high enough to 
allow the number density of the species in an excited state to be 
described by the Boltzmann distribution 



Q(T) 



(12) 



where and are the statistical weight and energy of the up- 
per state, respectively, Q{T) is the partition function, T^x is the 
excitation temperature, «x = «h,o^j:,o(1 + z)'^ with X^^ as the 
abundance extrapolated at the present time and the present den- 
sity of hydrogen atoms «h,o = Pcritf^B/'Wp ~ 2x10"^ cm"-'. We 
thus assume that the molecular abundances have reached their 
asymptotic limit at the redshifts of interest < 100 -200 on the 
large scales encompassed by the Odin beam. The possibility of 
a detection prior to a redshift of a maximum abundance is very 
limited due to the very low abundances and weak signals at even 
higher redshifts. On scales smaller than the Odin beam the abun- 
dances will certainly vary in regions with higher than average 
density or by the effects of the first stars. In reality, the popu- 
lation distribution also will depend on the molecular formation 
process, and the radiation field with possible population inver- 
sion (maser effects). 

The optical depth at the centre of the line can be calculated, 
assuming LTE and a Gaussian line profile, using 



i, Vl6;r3y3,Au Q(T) 



- l)d/,(13) 



where we have taken the stimulated e mission into account, and 
Vu/ is the frequency of the transition (cf. lPersson et alibOOTh . As 
customary we have converted the line width in frequency Av to 
a Doppler velocity width Av. The integration is performed over 
the path length of the cloud, L(z) -h - h- Note, that the path 
length determines the line width in the strictly linear phase (cf. 
Table IC.31 on-line material). 

4.2. Turn-around phase: 

During this phase the tendency of a perturbation to collapse un- 
der its gravity just balances its tendency to expand with the rest 
of the Universe. The perturbation will appear non-moving and 
all species from every part of the perturbation will emit from the 
same redshift. This will produce the strongest and the most nar- 
row resona nt lines with a line width determined by the thermal 
broadening dMaoli et al.Lll996l) 



Av 

V 



2 2ln2kTK 



7.16 X 10" 



(14) 



where m is the atomic mass and A is the mass in atomic mass 
units of the species. At low temperatures this implies a line width 
of a few km s"' at high redshifts for HeH* (cf. Table lC.3l on-line 
material). 

In addition to this, turbulence will also contribute with an 
unknown, and perhaps dominating, amplitude. A large region 
will also have a number of clumps which have narrow lines at 
slightly different velocities. The total line will be a superposition 
of these lines which will broaden the line. Finally, the cloud also 
consists of a number of smaller clouds that already have reached 
their turn-around which also may act as broadening agents. As a 
first order approximation in our analysis, we will therefore use 



a line width of 30 km s ' in our predictions of HeH^ resonant 
lines in the turn-around phase in Sect.|5] 

The integrated intensity of the resonant lines has not changed 
compared to the linear phase, and thus the amplitude, described 
by Eq. (|9]), will increase proportionally to the decrease in line 
width, which could be up to three orders of magnitude. 

As pointed out by iDubrovichI (f 19971) . a luminescence ef- 
fect may also produce lines when high energy photons are 
scattered and decay into several lower energy photons. The 
amplitude of lines produced by luminescence is given by 
Eq. ^ times an additional gain coefficient factor K. For sim- 
ple elastic scattering K=l a nd for luminescence K < 1000. 
iDubrovich & Lipovkal ( Il995h investigated luminescence pro- 
duced by excited rotational-vibrational (ro-vib) transitions of 
HiD^ which decay to the lowest vibrational state emitting via 
several rotational transitions. This effect is, however, limited to 
high temperature or high density regions. 

The turn-around and the beginning of the collapse phase 
are identifi ed as the most fav ourable for observations of reso- 
nant lines ( Maoli et"an Il996l) . It is therefore important to try 
to estimate at what redshift this will take place (full derivation 
in Appendix |B] on-line material). We assume that the power 
spectrum of the dark matter density fluctuations is of Harrison- 
Zeldovich type, which means 



O-M - 



AM 



2,1/2 



M 



-2/3 



(15) 



and that it can be normalised by the observed fluctuations 
a-M, = 10 '^ of the mass Ms =3.72 x IO'^Mq (cf. Eq. (15.13) 
iLongaiil I2OOI within the sound horizon at the last scattering 
surface where z lss = 1 090 giving 



(16) 



The turn-around redshift zta for a mass M which has 1 ctm over 
density, i.e. AM/M = ctm, can then be estimated according to 



1+Zta = 1.35(1 h-zlss) 



lO'^Me 
M 



2/3 



(17) 



Figure |4] shows the turn-around redshift z as a function of mass 
for one, three, and six ctm mass perturbations. Note, that the 
Harrison-Zeldovich power spect rum has the power ind ex n-1, 
and the latest WMAP results (Kom atsu et al.L l2009l) indicate 
that n =0.960^o'oi3- In our calculations we have taken the to- 
tal mass into account, but the baryonic to dark matter mass ra- 
tio is only 1/6. The star forming baryonic matter is predicted to 
subsequently fall into the dark matter potential wells caused by 
gravitational collapse. Note, that the collapse does not imme- 
diately follow when the object arrives at the turn-around point 
since the infalhng mass ha s difficulties to loose its energy (e.g. 
iDiemand & KuhlenLl2008h . 

To be able to quantify the signals from the proto-objects to 
observables we need to relate the turn-around mass to a linear 
size corresponding to an angular size of our beam. At the turn- 
around (1 + Ap/p)=(3;7r/4)^~5.55 and this should be inserted in 
Eq. (|7]i to find a mass-linear size relationship at the turn-around. 
We then use Eq. (I5]l - (l7]i and (fTTI l to find the relation between an- 
gular size vs. mass perturbations at turn-around redshift, shown 
in Fig. ID 
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Fig. 4. Turn-around redshift for one, three and six ctm mass over- 
dense primordial perturbations as determined by Eq. ( fTTT i. Note, 
that this describes the total mass, and the baryonic to dark mat- 
ter content is 1/6. Star and galaxy formation starts at later times 
since the baryonic matter does not collapse immediately at the 
turn-around redshift. 
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Fig. 5. The angular size vs. total mass at zta, assuming one, three 
and six ctm over densities. Note, that this describes the total 
mass. 



create a complex, multi-phase interstellar medium with a large 
range of densities and temperatures; up to 6 orders of magnitude 
at a given radius (Wise & Abel, 2008). 

The very energetic radiation from the first collapsed ob- 
iect s can for instance d issociate molecular hydrogen (e.g. 
Hai man et a ll lT997ll2000h . However, if the amount of H2 pro- 
duced in the cooling gas behind shock waves from the first 
supernova explosions exceeds the destroyed primordial H2 in- 
side the photodissociation regions, the first objects would have a 
net posi tive feed back on molecular production and galaxy for- 
mation. iFerrar i ([1998) concludes that multi-supernova explo- 
sions propagating through the interstellar medium can produce 
regions with a very high H2 abundance. For a wide range of 
physical parameters in such regions he found an H2 fraction of 
about 6x10"^. In addition, extremely energetic pair-instability 
supernovae (140-260 Mq) produced l arge amounts of dust at 
very early times (e.g. Schneider et al.L 120041) thereby allowing 
an increased molecular production including H2 and CO. All 
these effects influence the amplitudes of the primordial reso- 
nant lines and thermal emission and absorption. There are, how- 
ever, largejincertainties in the complex H2 and HD chemistry 
rilover & Abeii r2008). 

The resonant scattering will become suppressed when the 
density becomes higher than some critical density depending on 
species, since frequent collisions will cause thermal radiation to 
dominate over resonant scattering (Basu, 2007). The increasing 
numbers of collisions will drive the level populations and the 
excitation towards the lower (or higher) matter temperature, and 
should enhance the spectral line absorption (or emission). This 
effect will be greatest in molecules like H2 and HD, which have 
zero or small dipole moments and thus weak transitions. 



5. Sensitivity analysis 

Our observations have given us upper limits in terms of noise 
levels, which we now want to analyse using Eq. (|9]l. This 
equation is, however, expressed in intensity while sub-mm 
and radio antennas are usually calibrated according to the 
brightness temperature i n the Rayleigh-Jeans limit in terms of 
(lRohlfs& WilsonLl2004^ 



2k 



(18) 



Thus, the solution of the radiative transport equation with a 
constant source function, in our switched observations between 
signal and reference positions, becomes (cf. lPersson et al.L [20091) 



4.3. Non-linear collapse phase: 

When (1+ Ap/p) ^ (37r/4)^ the collapse will begin. The ampli- 
tude and line width now depend on the collapse rate in addition 
to the peculiar velocity. Depending on the ratio of the two veloc- 
ities resonant lines can appear in abs orption, emission o r have a 
double peak (a summary is found in lMaoli et all [19961) . During 
the initial stages the line widths are expected to increase com- 
pared to the narrow width during the turn-around. This is also the 
only evolutionary stage where the peculiar velocity is no longer 
required to produce resonant lines. 

The chemical abundances could also be substantially mod- 
ified during the collapse phase of a primordial cloud when the 
increasing temperature and density induce new chemistry. In 
such regions, molecul es will both be destroyed and produced 
(iPuv & Signorei Il996l) . The first Pop 111 objects will quickly 



ATI 



Th rjmb rjbi = [J(Tex) - J(Thg)] (I - e "") //mb Ihf 



(19) 



where Ar^ is the observed antenna temperature, ri,^b is the main 
beam efficiency, /^bf describes the beam filling (see Eq. (ISTT i be- 
low), Tex and Tbc, are the excitation and backgrou nd tempera- 
tures, and the radiation temperature, J(T) is CRohlfs & Wilsonl 
i2004l) 



J(T) 



hv 



1 



k 1 



(20) 



Note that the approximation is only valid if hv <K kT which 
is often used at frequencies < 100 GHz and temperatures much 
higher than the CMB. Figure IA.3I (on-line material) shows 
JiTcMB) as a function of frequency with T - 2.725 K. Around 
550 GHz, the approximation is no longer valid due to the low 
(CMB) temperature and high frequencies which lowers J(Tcmb) 
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Fig. 6. These plots are an example with the main goal to give an idea of the pattern of the primordial spectral lines: predicted 
antenna temperature Ar^ of resonant HeH^ optically thin lines within the 2'. \ Odin beam at z = 10, 30, 100, and 300. Turn-around 
phase is assumed for z = 10 and for higher z we assume linear phase with line widths as indicated in the legends. We also assume a 
low density hh = 2x10"^ and X(HeH^) = 10"'^. Note, that any increase of the density, abundance, or decrease of line width, 

directly increases the amplitude of the signals with an equal magnitude. If thermal absorption processes are present, they may be 
several orders of magnitudes stronger than the resonant lines. The amplitudes of the ro-vibrational lines are lower than 10"** mK 
until z ~ 300, when they show up around 200 to 300 GHz. They are, however, orders of magnitudes weaker than the rotational lines. 



to ~10"''. Accordingly, we use /(Tex) and not Tex in our calcu- 
lations. The 5cr noise level from Table [3] and |5] is used as Ar^ 

analysing the observations. 

The beam filling T]\,f can be estimated by (cf. iPersson et al.L 

120091) 



^2 _|_ 02 
^ beam 



(21) 



where 0^ is the FWHM angular size of the circular source, and 
^beam is the FWHM beam size. If the source is much smaller 
than the beam we will severely suffer from beam dilution by this 
factor We therefore assume that the primordial perturbation fills 
the large Odin beam (z/bf =1) thereby determining the minimum 
size of a perturbation to which we are sensitive (Figs. |4] and |5]l. 
In Table |C2] (on-line material) the physical sizes corresponding 
to the Odin beam as well as the step sizes in the 2006/07 obser- 
vations are given for a number of redshifts. According to Figs.|4] 
and|5] the Odin beam size is most sensitive to about 6x10'^- 



lO'"* Mq perturbations with turn-around redshifts of about 3-10 
(l-3cr). 

5.1. Resonant lines 

Equation (|9]l gives the intensity of the resonant lines. To ex- 
press this in Rayleigh-Jeans (RJ) brightness temperature we use 
Eq. ( fT9b and obtain 



n = J(Tcmb) (1 - e-') (3 - ay) cos ( 

c 



(22) 



Note, that the observed antenna temperature is lower by a factor 
m Vmh (Eq. (US). 

To estimate the maximum RJ brightness temperature pos- 
sible we assume that the lines in Eq. ( l22l l are optically thick 
(r > 1 ) since this is a limit after which the line will no longer 
increase its amplitude. The results are given in Table. |6] with 
an accuracy depending on the uncertainty of the assumed pe- 
culiar velocities. Obviously, even taking this uncertainty into 
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account, the resonant lines would be extremely difficult to de- 
tect at frequencies around 500 GHz and require noise levels or- 
ders of magnitudes lower than our observations. The Odin main 
beam temperature r^b-scale is related to the flux density by 
FylT^b = 2 600 Jy/K at 543 GHz. Figure |A3] (on-line material) 
shows J{Tcmb) X (3 - a), and demonstrates that we loose al- 
most two orders of magnitude in sensitivity around the observed 
Odin frequencies compared to frequencies below 200 GHz. This 
is a very limiting factor for detecting the already weak resonant 
lines, especially for the rotational-vibrational transitions which 
cannot be observed below 200 -400 GHz if they form around 
z = 200-400. The conclusion is therefore that the rotational- 
vibrational transitions of resonant lines most likely will be too 
weak for a detection even with future much more sensitive an- 
tennas. 

Table 6. Estimated largest possible brightness tem- 
perature Tb for optically thick resonant lines, and 
Vp(z)-(iQQ Hi -i-zkms"'. This velocity is typical for clusters 
while for a smaller mass it could be higher. Note, that these 
values are the same for all species and transitions. 



Obs. Freq. Maximum | I [mK] 



[GHz] 


z=10 


z = 30 


z= 100 


1 


1.6 


0.98 


0.54 


50 


1.5 


0.92 


0.51 


100 


1.3 


0.76 


0.42 


200 


0.64 


0.38 


0.21 


300 


0.24 


0.14 


0.078 


400 


0.071 


0.042 


0.024 


500 


0.019 


0.011 


0.006 


600 


0.005 


0.003 


0.002 



At frequencies < 100 -200 GHz the amplitudes of the res- 
onant lines are orders of magnitudes higher than around 500 - 
600 GHz, and hence a detection may be possible if the opacity is 
high at low frequencies. However, the resonant lines suffer from 
their dependence on the peculiar velocity which lowers their am- 
plitudes by approximately three orders of magnitude at z = 0, and 
even more at higher redshifts (Eq. (fTTT l and (l22li). Since thermal 
emission or absorption lines do not depend on peculiar veloci- 
ties, if existing, they have the potential to be much stronger than 
the resonant lines. 

Note, that these results also need to be corrected for beam 
filling, if the object is smaller than the beam, and for the beam 
efficiency to obtain the observable antenna temperature. For the 
Odin satellite this only lowers the amplitude by 10%, but for 
other antennas the reduction may be more than 50%. 

We also estimate the amplitude of optically thin resonant 
lines. As an example, with the main goal to show the pattern 
of the resonant lines as a function of redshift and to get a first 
rough estimate of the amplitudes, we consider the rotational 
and first branch of rotational-vibrational resonant HeH^ lines 
for four different redshifts, shown in Fig. |6] A single perturba- 
tion is assumed, with a size filling the 2'1 Odin beam. We have 
assumed a low density of hydrogen equal to the present one: 
nH = 2 X 10^^(1 + zy' at redshifts >30, and increased this density 
5.55 times at z = 10 (density at turn-around). We have used an 
HeH^ abundance relative to hydrogen of X(HeH^) = 10 and 
the background temperature follows the CMB at all redshifts. 
The adopted lines widths are narrow with AL' = 30kms ' at 
z = 10, and broad with Ay - 500 km s"' from z > 30 as indicated 



in the legends. This choice of densities and line widths arises 
from the fact that for z ^ 10, the large Odin beam corresponds 
to much larger physical sizes and masses than the turn-around 
objects (Fig. |4|, and at higher redshifts we are therefore most 
sensitive to objects in their linear phase. The peculiar velocity is 
described by Eq. ( fTTT ) with L'po~ 600 km s"' . 

Since the CMB temperature follows 2.725(1 -i-z)K, this im- 
plies that high redshift sources have an increased ability to pop- 
ulate higher molecular energy levels since the rising temper- 
atures give the photons more energy to excite the molecules. 
This is shown in Fig. |6] where high molecular transitions are 
increasingly populated towards higher redshifts. The rotational- 
vibrational transitions are not seen at all until z ~ 300, and are 
always much weaker than the rotational lines at all redshifts. 
It is also clearly seen that lines observed at frequencies around 
500 GHz always will be weaker than around 100 GHz; a clear 
disadvantage for our current observations. This could to some 
extent be com pensated by the p reviously mentioned lumines- 
cence process ( lDubrovichill997h but has not been considered 
here. 

We should here stress that ffie predictions of Fig. |6] are based 
upon rather conservative assumptions. Any increase of abun- 
dance and density, or decrease of line widffi, would directly 
increase the amplitudes of the HeH^ signals. Note also that 
ffie abundance is a function of redshift and that we here sim- 
ply have assumed that the abundance has reached a maximum 
value. The HeH^ abundance is indeed v ery uncertain with a pre- 
dicted abundance range of lO '^- 10-15 (iDubrovich etalilloosi: 
IVonlanthen eTall 12009). 

5.2. Thermal absorption lines 

The advantages of the Odin satellite are the absence of adverse 
effects from the terrestrial atmosphere and the tunable SSB re- 
ceivers, allowing a spectral coverage of a broad band. However, 
one severe disadvantage for detection of resonant lines with 
Odin is the high frequencies. To detect resonant scattering lines 
the best observing frequencies are below 100 -200 GHz consid- 
ering the rapid decrease of the radiation temperature as shown in 
Fig. lA.3l (on-line material). 

There are, however, other possibilities of detecting signals 
from the dark ages. Thermal emission or absorption depend on 
the competition between radiation and collisions as given by the 
radiative transfer equation 

'^^ " ~ir~ (e'"'..//*r„ _ 1 ~ ehv,,/kT,^ _ i) ~ ^ '^^ ' ^^^^ 

where T^ is the Rayleigh-Jeans brightness temperature, Tex is the 
excitation temperature of a transition and Tbg is ffie background 
radiation. Note that the latter two temperatures, as well as the 
opacity, have to be determined at the high redshift source, and 
that the frequencies on the right side in Eq.|23]are not the same. 
The excitation temperature depends on ffie radiation field, which 
is the CMB at high redshifts, and at lower redshifts it is domi- 
nated by the energetic radiation from the first star formation, and 
on the densities through collisions. The frequency in the first fac- 
tor is the redshifted rest frequency of the transition which is the 
same as the observed frequency, while the other factors include 
the rest frequency of the transition at the high redshift source. 

Equation (l23T l states that a strong signal requires a large 
temperature difference between the excitation and background 
temperatures, which can be produced by for instance collaps- 
ing primordial perturbations or the more rapid cooling of the gas 
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temperature compared to the CMB at even higher redshifts. If 
Tex ^ Tb„ Eq. (|23]) simplifies to 

/iVnhs 1 

The maximum Rayleigh- Jeans brightness temperature of an ab- 
sorption line can be estimated using Eq. (l24l l and assuming an 
optically thick line, which will cause the last factor to approach 
unity. The largest possible as a function of observation fre- 
quency when the background temperature is the CMB at the cor- 
responding redshift are given in TableQ We have also calculated 
the largest possible intensity when the excitation temperature is 
only 10% and 2% less than the CMB using Eq. (|23ll. Note, that to 
calculate the observed antenna temperature the brightness tem- 
perature needs to be corrected with the (unknown) beam-filling 
and beam-efliciency using Eq. ( fT9] l. 

Also note that if Tbg = Tcm&{z), the argument of the expo- 
nential functions in Eq. ( l24l i can be written as 

where Tq is the CMB temperature today. This implies that the 
maximum strength of an absorption line is redshift independent 
assuming an optically thick transition. In this case, at a chosen 
observation frequency the amplitude is the same for all species 
and transitions and is only determined by the known curve J{To) 
in Fig. |A3] (on-line material). 



Table 8. Estimated required background temperatures to pro- 
duce absorption lines with T\, = -lOmK at an observing fre- 
quency of 550 GHz assuming Tex <k T\,g, and an opacity of one 
or 10"^. Note, that since the observing frequencies are the same 
for all transitions the source redshift is different for every transi- 
tion. As a comparison, the CMB temperatures are listed for each 
redshift. 



Transition Redshift Opacity T\,g Tcmb 





z 


r 


[K] 


[K] 


n2(,z-u; 


lo.j 




JO 


J J 






10-3 


260 




H2(3-l) 


31.0 


1 


96 


87 






10-3 


430 




H2(4-2) 


43.4 


1 


132 


121 






10-3 


595 




H2(5-3) 


55.4 


1 


168 


154 






10-3 


756 




HD(4-3) 


18.1 


1 


57 


52 






10-3 


257 




HD(5-4) 


22.7 


1 


71 


65 






10-3 


317 




HeH+(4-3) 


13.4 


1 


43 


39 






10-3 


194 




HeH+(5~4) 


16.7 


1 


54 


48 






10-3 


240 





Table 7. Estimated largest possible absorption intensity T\, at dif- 
ferent observing frequencies for optically thick absorption lines 
assuming that the background temperature is the CMB. We have 
calculated this using Eq. ( l23l l and Tex - 0, Tex - 0.90 x Tcmb and 
Tex = 0.98 X TcMB- Note that these values are redshift indepen- 
dent. 



Obs. Freq. Max intensity | Tj, | 





rex = o 


Tex = 0.90 XTcMB 


Tex = 0.98 XTcMB 


[GHz] 


[mK] 


[mK] 


[mK] 


1 


2700 


270 


54 


50 


1700 


250 


51 


100 


1000 


210 


42 


200 


290 


97 


21 


300 


73 


33 


7.5 


400 


17 


9.1 


2.2 


500 


3.6 


2.2 


0.6 


600 


0.74 


0.5 


0.1 



put upper limits using the Odin observations depends on both 
the unknown background temperature and the redshift, in addi- 
tion to the opacity and excitation temperature. We give a very 
simple example of such an analysis in Table|8] We have here 
estimated the required background temperature (other than the 
CMB) in order to produce an absorption line with T\, = 10 mK 
at an observing frequency of 550 GHz for some molecular tran- 
sitions. We have used two different opacities and assumed that 
Tex T\,g (Eq. (l24l i). Such an absorption line would correspond 
to about a 5 cr detection for our 2009 deep search using 16 MHz 
channel spacing, or a 1 - 2 cr detection for our line surveys us- 
ing 16 MHz channel spacings (Table|5]l. We have also listed the 
CMB temperature at the corresponding redshift as a compari- 
son. A more sophisticated modelling, including excitation tem- 
peratures and estimated opacities, is demanded to produce more 
realistic results and to investigate how the intensities of different 
transitions vary with the physical conditions at different epochs. 



If the background temperature is not the CMB, the argument 
of the exponential function then becomes 

hVul /lVobs(l -Hz) 

As seen, the maximum absorption strength in this case depends 
on both the redshift and the unknown background temperature. 

The noise levels in our observations can be used to set up- 
per limits on possible absorption lines at the observed frequen- 
cies. The background source cannot, however, be the CMB at 
these frequencies since the maximum possible absorption inten- 
sity of such a signal is about 3 mK (Table|2l). This corresponds 
to about Icr of our best observations (deep searches), and even 
less for the spectral line surveys. Thus, the analysis in order to 



6. Summary 

In order to constrain cosmological models of star and struc- 
ture formation as well as the chemical evolution in the early 
Universe, we have performed spectral line surveys towards sev- 
eral positions without any known sources of emission in a search 
for primordial spectral lines from the Dark Ages. The first survey 
covered a broad band of 3 1 GHz between 547 - 578 GHz towards 
two positions with fixed reference-positions. The second sur- 
vey covered 1 1 GHz in the bands 542.0-547.5 GHz and 486.5 - 
492.0 GHz towards four positions. We also performed two deep 
searches towards one position with 543.100 and 543.250GHz 
as centre frequencies. No lines were detected, and thus the re- 
sults are upper limits in terms of noise level. Typical Icr values 
are 5-35mK in the 11 GHz survey, 14-90mK in the 31 GHz 
survey, and 2-7 mK in the deep searches over a 1 GHz band. 
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The maj or improvement made by the Odin observations 
compared to Ide Bernardis et alJ (1 19931) is the broad bandwidth 
covered allowing a wide range of redshifts to be explored for 
a number of atomic and molecular species. In addition, in the 
second survey we have taken into account the unknown sizes of 
the clouds by testing an observational strategy where we have 
observed towards four positions in a sequence with different an- 
gular distances between the reference and signal position. An 
important benefit of our observations is that we do no suffer from 
spectral line contamination from the terrestrial atmosphere. 

At low densities and in matter-radiation equilibrium condi- 
tions, the only expected signal is from resonant line scattering 
between CMB photons and matter moving with respect to the 
expansion of the Universe. These lines suffer, however, from the 
dependence of the low CMB temperature and the peculiar ve- 
locities of the moving primordial perturbations which is about 
10"^ at present and even lower at higher redshifts. In order to 
obtain an estimate of the highest possible intensities we assume 
optically thick resonant lines and no beam-dilution. The intensi- 
ties will then be at most of the order of a few mK at frequencies 
below 100 -200 GHz, and orders of magnitude lower at higher 
frequencies (5 fiK at 600 GHz). Since the lines most likely have 
T<^1 this will further lower their intensities with orders of mag- 
nitudes. 

If existing, thermal absorption lines on the other hand have 
the potential to be observable at both high and low frequencies. 
Such lines do not suffer from the dependence of peculiar veloci- 
ties of the moving primordial perturbations. The background ra- 
diation could also be considerably higher than the CMB, thereby 
producing stronger lines. The background could for instance be 
collapsing primordial perturbations, the first stars which are pre- 
dicted to form at z~20-30, or the remnants of pop 111 super- 
novae. These objects probably emitted large amounts of ener- 
getic radiation which foreground primordial perturbations could 
absorb. Also here we find the highest possible intensities by as- 
suming optically thick lines and no beam-dilution. Assuming 
that the CMB is the background radiation, we find that the largest 
possible intensity is 1 - 2 K below 100 GHz, and of the order of 
a few mK around 500 GHz. If the background radiation is higher 
than the CMB, the intensities will increase. Also in this case the 
opacity most likely is «: 1 which will lower the absorption line 
intensities, in addition to beam-dilution. 

The strength and line-width of the primordial lines depend 
on the evolution of the primordial perturbations. In the first 
(linear) phase they are very broad and weak, but become in- 
creasingly stronger and more narrow during the evolution of the 
cloud. The turn-around phase, which produces the strongest and 
most narrow lines, and the beginning of the collapse have previ- 
ously been identi fied as the most fav ourable evolutionary phases 
for observations dMaoli et al.Lll996l) . We have made a simple es- 
timation of the redshifts and sizes of the primordial perturbations 
at their turn-around phase, and find that the Odin beam size of 
2'.1 corresponds to a turn-around mass of about 6x10'^ Mq at 
a corresponding redshift of about 3 at a Icr level. Smaller per- 
turbations at higher redshift wiU thus suffer from beam-dilution. 
A beam-size of 40", which is approximately the beam-size of 
Herschel Space Observatory^, corresponds to a Icr turn-around 
mass and redshift of about 10'^ M© and z~ 10. 

The lowest rotational transitions of H2 will fall in the Odin 
band around 500 GHz from z~20-30, while the lowest transi- 
tions of other species such as HD, HD^, and HeH^ will fall in the 
band around 100 GHz from the same epoch. These transitions 
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may be searched for using ground based antennas - even though 
the foreground and atmospheric radiation will pose a problem. 

An important aspect of our work has been to test different 
observational strategies to prepare for our forthcoming obser- 
vations with the much more sensitive telescope and receivers 
aboard the Herschel Space Observatory launched on May 14, 
2009. The much lower noise level and the even broader band 
coverage with Herschel will increase the possibility of a de- 
tection. Other interesting facilities to consider in the searches 
for primordial molecules are for example the Atacama Large 
Millimeter Arra}0 (ALMA), the Combined Array for Research 
in Millimeter-Wave AstronomyQ (CARMA), the IRAM Plateau 
de Bure Interferometer and the IRAM 30-m telescope 0, and the 
Very Large Arra}|3 (VLA) at 20-45 GHz. Our observing meth- 
ods and resulting limits, paired with a sensitivity analysis taking 
into account the evolution of primordial perturbations, should be 
a valuable input to the planning of these observations. 

Spectral lines from primordial atoms and molecules may 
very well be the only way to probe the epoch of the cosmic 
Dark Ages and its end when the first stars formed. The search 
for these primordial signals is a very difficult task indeed, but a 
detection could be possible with the use of future facilities and 
would introduce an additional important way to discriminate be- 
tween models of the early universe as well as star and structure 
formation. 
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Fig. A.l. Spectra from weekend 1-5 vs. channels for Position 
A with B as reference from the 2006/07 survey in original shape 
and channel spacing 0.62 MHz. 



Channel 

Fig. A.2. Spectra from weekend 6-10 vs. channels for Position 
A with B as reference from the 2006/07 survey in original shape 
and channel spacing 0.62 MHz. 
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Fig.A.3. The blue thick line shows the radiation temperature 
J(To), and the red thin line shows JiTo) x (3 - a) as a function 
of frequency for Tq = 2.725 K. 
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Appendix B: The turn-around redshift as a function 
of density fluctuations on tlie last scattering 
surface 

Let the unperturbed density at the last scattering surface (LSS) 
be plss = PniLss + PAlss + PiadLss' whcre p^^^j.^ is the matter 
density, Pa^ss is the dark energy density and Pradtss '^he radi- 
ation energy density. Since the matter density is very domi- 
nant at LSS Plss ~ PniLss- This implies that the critical density 

PcritLss =PmLss +PAlss +PradLss + Pki^s ~ PniLss sinCC WC aSSUmC 

that the curvature density pkLss is zero. The density parameter at 
the LSS then becomes Qlss = PLSs/PcritLss = 1- 

Let a perturbed density at LSS p^gg with the (large) scale 
flLSS be defined by 



Plss = Plss + ApLss , 



(B.l) 



where ApLss is the density perturbation. If ApLss is po sitive, the 
turn-around scale a'^^^ of the perturbed region is (cf. iLongaiii 
[2008. Eg. (16.2)) 



Ota — '^LSS 



where 



LSS 



„, _ Plss + ApLss _ , ApLss 
'"'lss ~ - i + 



Plss 



Plss 



(B.2) 



(B.3) 



If ApLss is negative, there is of course no turn-around. 

The scale oja in the unperturbed region at the turn-around 
time t - fxA is then (cf. Longair (cf. Longair, 2008, Eq. (16.3)) 



ajA = 5.55'/^ ■ fl^A = 5.55'/^^ ■ alssttt 



1/3 , 



^■'lss 



"lss-1 



(B.4) 



Since ota = (1 + Zta) ' and alss = (1 + Zlss) ' we have 



sound horizon at the last scattering surface where zlss - 1 090. 
Equation (IB.8I 1 is then normalised to 



1 0-4/^^1 



2/3 



(B.IO) 



Expressing AM/M in Eq. (IB. 7b in terms of ^ ■ <tm, where ^ is 
a sigma measure, we get 



(l+ZTA) = (l+ZLSs)-5.55-'/^-^-crM 

Note that ^ > 0. 

Inserting Eq. (IB.lOl l and Ms we get. 



1 +ZTA = (1 +zlss) ■ L35 - ^^ 



M 



(B.ll) 



(B.12) 



This derivation is valid for dark matter Baryon structure for- 
mation, also driven by gr avitat i onal f orces, has as well been 
affected by Silk damping JSilH 1 19681) and sound waves. This 
complication introduces a minor modification to the power spec- 
trum of the density fluctuations (page 412 in Longair, 2008), but 
approximatively our derivation should hold also for the bary- 
onic matter. More worrisome for the Harrison-Zeldovich sce- 
nario are perhaps the findings of massive galaxy clusters at 
high redshifts, for instance XMMU J2235.3-2557 at z = 1.4 
with a estimated projected mass of the cluster within 1 Mpc 
of 8 .5 ± 1 .7 X lO'-^Mo (Je e et aU 120091) which according to 
Eq. ( IB.I2I 1 constitutes a 14cr event. 

It should perhaps be noted that masses on the order of 10^ - 
10^ Mq have turnaround redshifts on the order of lO"' - 10"*. This 
means that black holes of this size could form shortly thereafter 
if the collapse is not halted by induced angular momentum. 

Appendix C: Tables 



(1+Zta) = (1+zlss)-5.55 
If ApLss <^ Plss we have 
(1+Zta) = (1+zlss)-5.55 



-l/3_ 



ApLSS 



Plss + ApLss 



(B.5) Table C.2. Physical sizes corresponding to the Odin beam and 
the different angular steps in the 2006/07 observations. 



-1/3 APLSS 
Plss 



the mass within the scale a, so that 



(1+Zta) = (1+zlss)-5.55 



_1^3 AMlss 

Mlss 





Angular 






Redshift z 








(B.6) 


size 


z = 5 


10 


25 


50 


100 


200 


300 


where M is 


en 






Physical 


size [kpc] 








2:1 


820 


540 


270 


150 


79 


41 


28 




5 


2000 


1300 


640 


350 


190 


98 


66 




10 


3 900 


2 600 


1300 


700 


380 


200 


130 


(B.7) 
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3 900 


1900 
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560 


290 


200 
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12000 


7 700 


3 800 


2 100 1 100 


590 


400 



Conside r now density fluctuations of the Harrison-Zeldovich 
type (cf. lLongaiitl2008l page 390-392), where 



0"M 



(B.8) 



Note that this relation may be expressed as 



(B.9) 



assuring that black holes are not formed excessively on small or 
large scales. 

Assume now that Eq. ( IB. 81 ) can be normalised by 
the observed fluctuations ctms - 10 of the mass 
Ms = 3.72 X 10'^ (cf. rLongaiiil2008l Eq. (15.13)) within the 
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Table C.l. The ranges in redshifts covered by our observations for the five lowest rotational lines of H2, HD, and HeH^. 



Species 


Frequency 


Transitions 


E„ 


A-coefF 




AC [GHz] 


AOS [GHz] 


AOS [GHz] 


AC [GHz] 




[GHz] 


/ 

•J LI 




[K] 


[s-'] 




486.5-492.0 


542.0-547.5 


547.0 


-563.0 


563.0-578.0 




10 621 


2 


~0 


510 


2.94x10" 


11 




1 8 60 - 1 8 40 


18.42 


-17.87 


17 87- 17 38 




17 594 


3 


-1 


1 015 


4.76x10" 


10 




^1 4(S-31 14 


31.16 


-30.25 


^0 7'S _?9 44 




24410 


4 


-2 


1 681 


2.75x10" 


-9 


2 = 49.18-48.61 


44.04-43.58 


43.63 


-42.36 


42.36-41.23 




31011 


5 


-3 


2 503 


9.83x10- 


-9 


z = 62.74 -62.03 


56.22-55.64 


55.69 


-54.08 


54.08-52.65 




37 348 


6 


-4 


3 474 


2.64x10- 


-8 


2 = 75.77-74.91 


67.91-67.22 


67.28 


-65.34 


65.34-63.62 




Z D / J 


1 


-0 


ITS 
IZo 


z 1 n- 
J.jZX iU 


-8 


■7 — A A 'XI 


Q Q/1 a on 


3.89 


-3.75 


J. IJ — J.KiJ 






2 


- 1 




c nc w 1 n- 
j.UjX lU 


-7 




o.o4 — 0. /4 


8.75 


-8.47 


AH Q 07 

0.4 / — o.ZZ 




7 952 


3 


-2 


766 


1.80x10- 


-6 


2=15.35-15.16 


13.67-13.52 


13.54 


-13.12 


13.12-12.76 




10518 


4 


-3 


1271 


4.31x10- 


-6 


2 = 20.62-20.38 


18.41-18.21 


18.23 


-17.68 


17.68-17.20 




13015 


5 


-4 


1895 


8.35x10- 


-6 


2 = 25.75-25.45 


23.01-22.77 


22.79 


-22.12 


22.12-21.52 


HeH+ 


2010 


1 


-0 


96 


0.109 




2 = 3.13-3.09 


2.71-2.67 


2.67 


-2.57 


2.57-2.48 




4009 


2 


-1 


289 


1.04 




2 = 7.24-7.15 


6.40-6.32 


6.33 


-6.12 


6.12-5.94 




5 984 


3 


-2 


576 


3.75 




2=11.30-11.16 


10.04-9.93 


9.94 


-9.63 


9.63-9.35 




7 925 


4 


-3 


956 


9.14 




2=15.29-15.11 


13.62-13.48 


13.49 


-13.08 


13.08-12.71 




9 821 


5 


-4 


1428 


18.1 




2=19.19-18.96 


17.12-16.94 


16.95 


-16.44 


16.44-15.99 



Table C.3. Line widths in the linear phase for a perturbation corresponding to the Odin beam at a number of redshifts. 



Redshift Temp Size" Mass" Fract. width Line width* Width at 500 GHz 



2 T L M Av/y AwLinear Av 

[K] [kpc] [Mo] [kms-'] [GHz]] 



10 


30 


540 


4xl0'2 


2x10-^^ 


720 


1.3 


20 


57 


320 


6xl0'2 


4x10"' 


1 100 


2.1 


50 


139 


150 


9xl0'2 


7x10-3 


2 000 


3.7 


75 


207 


100 


9xl0'2 


8x10-3 


2 500 


4.5 


100 


275 


79 


lxlO'3 


1x10-2 


2 900 


5.3 


150 


411 


54 


1x10" 


1x10-2 


3 600 


6.6 


200 


547 


41 


lxlO'3 


1x10-2 


4 200 


7.7 



"For an angular size of the Odin beam, 2'. 1 calculated with Eqs. (|5) - l|7]l. *The line width during the linear evolution, Eq. l[8j. 



